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Abstract 

Collective plasma excitations of optically dressed Dirac electrons in single and double 
graphene layers are calculated in the RPA approximation. The presence of circularly polarized 
light gives rise to an energy gap Eg between the conduction and valence energy bands. Its value 
may be adjusted by varying the frequency and intensity of the light, and may reach values 
of the gap reported for epitaxially grown graphene and far exceeding that caused by spin- 
orbit coupling. We report plasmon dispersion relations for various energy gaps and separation 
between graphene layers. For a single graphene sheet, we find that plasmon modes may be 
excited for larger wave vector and frequency when subjected to light. For double layers, we 
obtained an optical and phonon-like mode and found that the optical mode is not as sensitive 
as the phonon-like mode in the long wavelength limit when the layer separation is varied, for a 
chosen Eg. The dressed electron plasma although massive still has Dirac origin giving rise to 
anomalous plasmon behavior upon crossing the ft) = hvpq cone. 
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I. Introduction 



Recently, several groups have addressed the effects which a circularly polarized electromagnetic 
field (CPEF) ,12 spin-orbit interaction (SOI) 1^ in suspended graphene or the sub-lattice symmetry 
breaking (SSB) by an underlying polar substrate!^ in epitaxial graphene may have on the energy 
band structure of a graphene sheet. In all cases, a gap is opened between the valence and conduc- 
tion bands as well as between the intra-band and inter-band electron-hole excitation continuum of 
the otherwise semi-metal Dirac system .1^^ As a matter of fact, the interplay between the single- 
particle excitations in the long wavelength limit results in dielectric screening of the Coulomb 
interaction which produces an undamped plasmon mode that appears in the gap separating the two 
types of electron-hole modes forming a continuum. By this we mean that a self-sustained collective 
plasma mode is not supported by exciting either the intra-band or inter-band single-particle modes 
only. Furthermore, although the Dirac electrons near the K-points now acquire a non-zero effec- 
tive mass, one still cannot produce a long wavelength plasmon mode as in the two-dimensional 
electron gas (2DEG) by intra-band excitations for intrinsic graphene. 

Plasmon excitations in semi-metallic graphene were studied extensively theoretically in single 
and multilayer graphene. However, their experimental observations are limited so far to either in- 
elastic scattering of highly energetic particles (faster than 10^ m/s)ox indirect observation through 
plasmoritons (a quasiparticle formed by resonant coupling of renormalized hole states with the 
plasmons). Plasmons in graphene gapped by SOI or by a polar substrate have been investigated 
by Wang and Chakraborty ^ as well as Pyatkovskiy ^ Because the SOI is small (~ 0.05 meV), 
its effect on dielectric screening of the Coulomb interaction requires either high resolution spec- 
troscopy or low electron density in the conduction band to be observed. The effect of SOI on the 
plasmon excitations in a 2DEG were also studied by several authors On the other hand, the 
gap between the valence and conduction bands near the Fermi level when circularly polarized light 
is applied to graphene could be manipulated, unlike the intrinsic gap due to SOI, depending on the 
intensity and frequency of the light. For a CO2 laser with a 70 W output power and wavelength of 
10 ji, the energy gap is about 100 meV, compared with 0.08 meV induced by SOI .12 The exact gap 
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opening mechanism in epitaxial graphene is still under debate The problem is rather technical 
and arises from the fact that the graphene sample sits on top of a buffer layer which provides an 
additional mid gap level, thus obscuring the exact energy dispersion curve and requires numerical 
ab initio calculations This ambiguity stimulated discussion regarding the symmetry breaking 
gap [H versus the effects due to electron-electron interaction .ISlThis also not only indicates the im- 
portance of theoretical investigation of plasmon in gapped graphene, but also an alternative means 
of providing the gap. Since SOI is too weak to observe its effect, dressing of the Dirac electrons 
with CPEF is seems to be an alternative worthy of further investigation. 

The value of this dynamic energy gap is comparable with that induced by a polar substrate, but 
contrary to the latter, it is highly adjustable. As a matter of fact, those two gap-opening mechanisms 
may even work conjunctively. In that case, the circular polarization of the electromagnetic field 
may no longer be required provided that it is tuned to be in resonance with the SSB gap. A possible 
experiment to measure the plasmon frequencies in single and double layer graphene could employ 
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electron energy loss spectroscopy (EELS) ,^^211121 as shown schematically in Fig. [figure][l][]l 



sentially, this would be a pump-probe, with the electron beam being the probe effected by the 
polarized light whose role is to open the gap. The Raman shift of the scattered electrons provides 
both particle-hole and plasmon excitation frequencies, which are usually characterized by their 
spectral weight, a quantity that depends on the transferred energy hco and momentum hq. 

The remainder of this paper is organized as follows. In Sec. II, we present the model Hamil- 
tonian for dressed Dirac electrons by a CPEF as well as corresponding energy eigenvalues and 
wave-functions. These results are then used for calculating the polarization function and plasmon 
dispersion relation in the RPA. Our numerical results are presented in Sec. Ill but, for convenience, 
we briefly summarize the highlights of our calculations. For a single graphene layer, in the presence 
of CPEF, we calculate analytically the real and imaginary parts of the noninteracting polarization 
function , whose imaginary part provides the particle-hole spectral weight. As a piecewise con- 
tinuous function, it is expressed in terms of region functions which provide a convenient way for 
determining the Landau damping of the plasmon excitations. By increasing the intensity of CPEF, 
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the Dirac electrons become more dressed by photons, providing them with a larger and larger ef- 
fective mass. Regions in £0 — ^ space, which are non-existent for conventional Dirac electrons, 
emerge thereby making the plasmon behavior to resemble that of conventional 2DEG. However, 
the resemblance is limited, and we show there is no complete crossover to 2DEG behavior. The be- 
havior of the plasmons is definitely determined by the presence of the Dirac cone CO = Vfq, where 
vp is the Fermi velocity. This indicates that although the dressed electron plasma is massive, it still 
has Dirac origin. We also investigate graphene double layers subjected to a CPEF with two plas- 
mon branches given as resonances of the real part of the inverse dielectric function. Two cases are 
investigated: the symmetric configuration when both layers are gapped and the asymmetric case 
when the gap is induced on one of the layers. The latter case may be obtained by compensating 
CPEF with the SSB on the given layer. The main derived effect on the acoustic plasmon branch 
is that this mode shows a gain for the symmetric case, but a loss for the asymmetric case. This 
applies for a sufficiently intense laser field, over a restricted range of frequency and wavelength ^ 
In Sec. IV, we offer some remarks concerning experimental observation of the plasma excitations. 

II. Model and Formulation of the Problem 

Let us consider several graphene layers separated by distance d. One of the layers may be epi- 
taxially grown on a SiC substrate. The structure is optically pumped by an idealized single mode 
CPEF characterized by an associated energy Nohcoo where A^o is the average number of photons 
in the optical field. It is assumed to be a function of the external pump intensity. We also take it 
that all graphene layers are in the node of the optical field and we can neglect retardation effects. 
This is valid provided that d <^ Itzc/coq. The objective is to find the plasmon dispersion branches 
(0±{q) for double layer graphene. Those are given by the poles of the inverse dielectric function. A 
single layer may be described by the same formalism by setting the inter-layer distance to infinity. 
By setting the origin z = at one of the layers, we may follow the procedure presented in Ref.l^ 
to obtain the inverse dielectric function as: 
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e-\z,z';q,0})^5{z-z')+ £ 5{z-U-mvc{jdj'd)U^>{q,o}) 



(1) 



where 5(z — z') is the Dirac delta function and polarization matrix is 



^ nil ni2 ^ 



V 



(0) \ 



It is expressed in terms of the generaUzed dielectric function 



(2) 



e{q,CO) = l-Vc{q)\^ 



n(o) , n^o) 



+ 



[l-e^^'^')vM<nf^]} . (3) 



where Vc{q) = 27ue^/ Esq with = 4;reoe/„ £h being the average background dielectric constant. 
The plasmon resonances correspond to the poles of the inverse dielectric function which are given 
by the solutions of e{q, (d) = 0. Additionally, the plasmon frequency incorporates its a relaxation 
rate so that CO ^ co + irf. Also, Oyy^ (q^co + irj) is the noninteracting polarization function for the 
7* graphene layer. The inter-layer polarization function nj2^ is not included in our formalism due 
to negligible interlayer tunneling. 

The key ingredient of the above equations is the polarization function. This requires specifics 
from the model given by the Hamiltonian at the two inequivalent K, K' points in the graphene 
dispersion: 



H = J^jc + J^D (4) 
Jtjc = h(Ooa^a+ ^^(73 - ^= ((7+fl + (7_fl-f) (5) 

Mb = hvF ii<y+ + 0-)kx±i (a_ - C7+) ky) . (6) 



Here, the Jaynes-Cummings Hamiltonian (??) is described by the interaction third term governed 
by the amplitude of the electric filed Eq = ^JNoHcoq/ £qV , and V is the quantization volume of 
the CPEF. The energy of an electron in rotational motion induced by the field is denoted as Wq = 



5 



IvfcEQ/coQ ~ v^A^ and assumed to be much smaller than the energy of the optical filed itself, i.e., 
Wq/Tkoq <^ 1. The second term in ?? takes account of the breaking of the inversion symmetry 
between A and B sub-lattices. Also, A/i^^s stands for the difference in on-site energies for the 
sub-lattices. It may be induced if the graphene is grown epitaxially on a substrate, and may acquire 
values as large as 260 meV (4775 nm), see, for instance, Zhou et aL^^and references therein. On 
the smaller scale (~ 10"^ meV) the symmetry breaking term accounts for spin-orbit interaction. 

The Dirac part of the Hamiltonian is given by ??, with ± standing for K, K' points correspond- 
ingly. The operators involved in the above Hamiltonian are acting in the joined electro/photon 
space via the following relations: 



a+ = \t){i\; a- = |i)(t|; 

a+|t,A^o)=0; a+| i,A^o) = I t,A^o) 

a_|i,^o)=0; a_|t,^o) = |i,^o) 

a\ It.No) = v^l it,^o - 1) ; ^'^l it,A^o) = VN^\ it,A^o+ 1) 

with III) denoting the Dirac pseudo-spin basis. The Jaynes-Cummings part can be readily diago- 
nalized by switching to the dressed electron basis 



^ |+,A^o) ^ 



v 



|-,A^o) 



-sin(4>^o/2) cos {^N J 2) 



\ 
I 



cos(4>^„/2) 



sin(4>w„/2) 



2Q. 



IQ. 



'No 



2Q. 



'No 



lt,A^o) 

11,^0+1) 



(7) 



where = + '^Wq{No + l)/No and the detuning between the microcavity and the symmetry 
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breaking induced gap is denoted by A = hco — A/i^ g. In the new basis, the Hamiltonian ?? within 
the small dressing regime Wq <^ A assumes the following form 

H = Woh(Oo+^03+^^D (8) 

where Eg = ^Jw^~+~^ — hcoo is the energy gap between valence and conduction bands. The 
gap governs the nature of dressed Dirac electrons metal-insulator transition.'^' One may obtain the 
solution of the above equation in a straightforward way, giving 



v 



( 



1 



hvfk 



\ 



(9) 
(10) 

(11) 



It is now a simple matter to show that conventional Dirac fermions correspond to the off- 
resonance limit of the dressed states. Given the dressed states, the noninteracting polarization 
assumes the Lindhard form 



TlJ]{q,(0)=Tlf\q,(0)+nf^{q,0i) (12) 

U^]^\q,(0)=-'- Z jkdkJ^-{Kq,m); U^^^\q,(0) = )^ E kdkJ^^P {k,q,CO) 

a=±o a,j8=± 

J^'^Plka (0)-^Ud • g ^^'4k-(k+q)+(£,/2)^ \ 



The two terms in ?? stand for the inter- and intra-band polarizations, respectively. Also, /i is the 
chemical potential whose zero value is in the middle of the gap (Eg). In calculating the polariza- 
tion, we assumed that for chosen A^o^ each point on the dispersion curve might be occupied 



7 



by no more than one electron. This assumption is valid for small energy gap. Otherwise, the 
quasi-equilibrium density matrix acquires non vanishing off-diagonal elements thus substantially 
complicating the statistical properties of the dressed states. The noninteracting polarization ?? in 
gaped graphene has already been calculated on both the real and imaginary frequency axes'^'^ 
which essentially yield the same result. Since it plays a crucial role in this paper, we give its full 
form along the real frequency axis as 



Uf\q, £0) = + , X [ (13) 



(/G>(.^i,_)-/G>(^i,+)))l< + (G<(.^i,_)+/G>(xi,+)))2< + 
(G<(xi,+) + G<(xi,-))) 3< + {G<{xu-) - G<(xi,+)))4< + (G>(xi,+) - G>(xi,-))) 1> + 
(G>(xi,+) +/G<(xi,-))) 2> + (G>(xi,+) - G>(-xi,-) - i7t{2-xl))) 3> + 
(G>(-xi,_) + G>(.^i,+) - i7ti2-xl))) 4> + (Go(.^i,+) - Go(xi,-))) 5>] 

Here, the following notations for the region functions have been introduced: 



xq = J 1 + 



G<{x) 
G>{x) 
Go{x) 



p2 
^8 



hvfq 

h\l{q±kFY + {Ej2y 



xy x^ — x^ — (2 — Xq) cos~^{x/xo) 



x-i/x^ — Xq — (2 — Xo)cosh ^{x/xq), 



x-y/jc^ — Xq — (2 — Xo)sinh ^{x/J—Xq 
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and the regions are shown in Fig . | [figure] [2] [] 2| and are defined as 



1< = (iJ. - X2- - hco) 

2< = G {-hco - n + X2~) 6 {hcj + n - X2-) 6 + X2,+ - hcj) 

3< = 6{-iJ.+X2--fi(o) 

4< = {h(0 + iJ.-X2,+)6 {hvfq-h(o) 

1> = {2kf - q) 6 {hco - X3) 6 (iJ. + X2- - hco) 

2> = 6 {hco - IJ. - X2-) {iJ. + X2^+ - hcj) 

3> = 0(^CQ-/i-X2,+) 

4> = e{q-2kf)d{h(O-X3)0{il+X2--h(o) 
5> = d{hco-fiVFq)9{x3-hco) 

For the purposes of optical spectroscopy, it is advisable to study the long wavelength limit by 
expanding the polarization function in ?? in powers of q. For small q/kf, the undamped plasmons 
exist in the spectral regions 1> and 5>. By expanding the polarization up to the fourth order in q, 
we obtain 



2fi—hco 
2fi+hco 



Ahjico + [Ej + h^O}^) log 



2fl—hC0 



For small frequency ft), the following identity applies, (fi^ + ^^ft)^) log 
which we employed to further simplify the polarization (??) as 



2^+hco 



1 



>5 5>,4<; 

1<,2<. 

(14) 
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(l< + 2<) 



Note that smallness of (O is automatically assured in the regions 4<, 1<,2<. In these regions, 
the dominant excitations are of the particle-hole type, due to nonzero imaginary part of the po- 
larization. In the long wavelength limit, the plasmon frequency for a single graphene layer in the 
vacuum is obtained by solving for the zeros of e{q, (O) in ??, along with uf} = {q, (o) given 
by ?? and 1122'' set to zero, yielding 



(ol = q^{Eg,(Op) , (16) 



where we introduced the plasmon region function 



(17) 

In the above expression we used /Apieohvp — 2.5 It is also straightforward to show that the 
region function disappears in the region 1 <, thus supporting Landau damped plasmons in the region 
2< only. We also note that for conventional graphene and the 2DEG, the plasmon region function 
?? is replaced by constants Ifie^/h^Ss and m*Vpe^/h^£s, respectively. 

For the symmetric double layer configuration, one obtains two undamped and two Landau 
damped plasmon branches whose frequencies in the long wavelength Umit obey 

col = ^e"^'^ (e«^ ± 1) j3^{Eg, (0±) (18) 
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where ± stands for the bonding/anti-bonding plasmon modes, when the polarizations in the two 
layers work together/against each other in producing the charge density oscillations. Two inde- 
pendent parameters are qd and Eg. That is the interplay between the inter-layer distance and the 
plasmon region function determines the plasmon branches. In the limit of large inter-layer separa- 
tion, qd^ I, both branches are reduced to ??. Conventionally ,'22 those regimes are referred to as 
intermediate and weak coupling limits. We avoid using these terms so as not to confuse with the 
light-matter coupling. On the other hand, for small inter-layer separations, ^ 1, we have 



(ol = q0^{Eg,G)±){2-qd) (19) 
(ol=q^^{Eg,(0±)d (20) 

Due to linear dependence of the undamped plasmon frequency on the wave vector q, the anti- 
bonding mode is often referred to as acoustical branch. The other undamped plasmon branch 
has larger frequency and is named the optical branch. To the leading order in q, this branch is 
independent of inter-layer separation and is completely determined by the plasmon region function. 
It is always higher in frequency than the acoustical branch and usually (but not always) experiences 
Landau damping in the long wavelength limit. 

The asymmetric case is also easily obtained in the long wavelength limit, i.e., 

col = ^^e-'^^ (^e'?^(^i,± + ^2,±) ± V^e2^^(^i,±-^2,±)2+4^i,±^2,±) , (21) 

where we have introduced the shorthand notation ^i(2),± = '^(^^,1(2)5 ®±)' with £^gj(2) denoting 
the gap induced on the first (second) graphene layer. In the small inter-layer separation limit, these 
eigenfrequencies become 
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0)1= q (^2+ + _ M - 1)) ('^1,+ + (22) 

£0^ = 2J^2^i,_^2,- {^\,- + <^2,-)"^ ■ (23) 

When one moves away from the long wavelength limit, e{q, ft)) = from ?? in conjunction with 
(??) become nonlinear, with no closed-form analytic solutions. Consequently, numerical simula- 
tions are needed. Our numerical results are the subject of the next section. 



III. Numerical results and Discussion 

In our numerical calculations, we scaled energies in units of the chemical potential /i and the 
wave number in units of the Fermi wave number kf. The frequency of the electromagnetic field 
along with its intensity determines the energy gap Eg between the valence and conduction bands. 
Consequently, this is the only parameter we use to specify the role played by the external laser 
pump field. All our calculations for the collective plasma excitations were carried out at T = ^. 

We begin this section by considering the plasmon dispersion for conventional Dirac electrons 
in graphene obtained by solving e{q, ftj) = with the use of ??. Our results are presented in 



[figure] [3] []3 ^a.l) and are superimposed onto the particle-hole modes by plotting 3m Yl^^\q,a)). 
This is done in order to identify the regions where there is Landau damping. Clearly, the zeros of 
?? yield both damped and undamped plasmon excitations. In order to isolate the undamped modes, 
it is convenient to consider the plasmon spectral weight —(1/ 7r)3m Tl{q^ (O — ijq), where Tl{q, ft) — 
ijq) = n^*^) {q, ft) — ijq) /£{q,(0 — ijq) . The plasmon relaxation rate is determined self-consistently, 
starting with the phenomenological rate whose value hereafter is chosen as fiy/ ji = 10^^. If not for 
the finite relaxation rate, the plasmon spectral weight is nonzero only in the particle-hole regions ^ 
Below, we show that this condition applies only for the single layer configuration. The concept 
is demonstrated in [figure] [3][]3| ^b.l). The naturally damped plasmon branch in the 1> triangular 



region becomes strongly Landau damped once it enters the particle-hole mode region 2>. 
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By dressing Dirac electrons with circularly polarized photons, they acquire an effective mass 
and a gap Eg is opened between the valence and conduction bands. The gap is determined by 
the intensity and frequency of the laser field. The presence of Eg in turn induces a gap in the 
particle-hole mode region for q/kp ~ 1. The particle-hole free region 5> widens with increasing 
Eg. Its contribution to the plasmon region function also grows, as indicated by the second term 
in ??. Meanwhile, the Dirac contribution of the plasmon region function is decreased through the 
first term in the same equation. Schematically, the opening of the plasmon window by the electron 



dressing is shown in [figure] [2] []2 There is one more region in the particle-hole continuum which 
gets opened by the energy gap Eg, namely region 4<. We give it special attention since this region 
provides the main difference between dressed Dirac electrons and conventional graphene. In the 
latter case, this region simply merges with 5>, since the (O = fiVfq line has no physical meaning. 
However, in gapped graphene, this region is very special due to the possible natural plasmon 
gain existing in it. Indeed, according to ??, the polarization acquires an imaginary part provided 
analytical continuation CO — t- £0 — zy from the real axis into the lower complex half-plane 7 > 0. 
The imaginary part of the polarization becomes a weighted delta function, centered at the plasma 
frequency 3mn(^, (O) = 3m(l/(vc(<?)e(<?, COp))) = ^gn{Zm£{q,(0p)){n/vc{q))5{^e£{q,0}p)) = 
nW {q,COp)d{(0- (Op). 

In the region 5> the plasmon weight is the conventional positive one given by W{q, (Op) = Wo5> 
, with Wo = \/v^{q)-^'^eTl^^\q,(o)\(o^co^^) > 0. This behaviour can be seen on the lower insets of 



[figure] [3] []3 (a.3). However, when the plasmon branch enters the region 4> the plasmon weight 



may flip its sign W{q,(Op) = — Wo4>, since 3Tl^^\q,(0p) < as indicated on the upper inset of 



[figure] [3] []3 ;a.3). Interestingly enough the change in sign of the plasmon weight requires the 



plasmon branch to cross ft) = hvfq line. On one hand, the gapeless graphene does not show the 
crossover from the plasmon loss to the gain since the above condition is never satisfied. On the 
other hand, if one is to increase energy gap (Eg/ /i ^ 2) so that the plasmon branch is completely 
forced into the 4< region, then the plasmons also show only the loss and the positive plasmon 
weight W{q, (Op) = WqA^. For even larger energy gap, the plasmons are forced into region 2< and 
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become fully Landau damped. The plasmon dispersion for a single layer of gapped graphene is 
presented in |[figure] [3][]3[ a.2), (a.3). Numerical results shown in |[figure][3][]3[ b.2), (b.3) demon- 
strate the predominant plasmon gain for short wavelengths. In the long wave length limit the 
plasmon branches are always well behaved. Therefore, at this level of theory, it is impossible to 
say if the anomalous plasmon behaviors has some physical meaning or being just an indicator of 
RPA failure at larger q. 



Turning now to the double layer configuration, we present in [figure] [4] []4 through [figure] [7] []7 



the plasmon dispersion and damping for the symmetric and asymmetric cases. In these plots, the 
panels (a.l, a.2, a.3, a.4) correspond to zero energy gap for chosen inter-layer separation. In the 
presence of an energy gap, we note that one of the plasmon modes (the acoustic branch (o ) may 
be forced into region 4< for values of Eg/ ji as small as one (See ??(b.l, c.l)). As far as the sym- 
metric configuration is concerned, this mode demonstrates gain, while the asymmetric case shows 
losses in that region as may be deduced from its spectral weight in ?? (b.l, c.l). 

For the symmetric double layer, there are two Landau damped plasmon modes in the region 
2< . One of them runs just along 2< , 4< interface. While the other starts along 2< , 1 < interface and 



then eventually merges with the first branch at short wavelength, as demonstrated in [figure] [4] []4 



The frequencies of these plasmon branches are almost independent of the inter-layer separation 
d. The asymmetric configuration in effect causes one of these damped modes to become purely 



acoustic with frequency CO ~ qvf over a wide range of wave vector q, as shown in [figure] [5] []5 



Let us now confine our attention to the plasmon branches whose frequencies in the long wave- 
length limit coincide with those of the bonding ((0+) and antibonding (CO-) modes. For large 
inter-layer distance satisfying kfd ^ 1, in the absence of an energy gap. Eg = 0, these modes con- 
verge onto the line CO + Vfq = 2ii/h. By increasing the gap for the symmetric double layer, there 
emerge Landau damped plasmon branches in regions 2>,4> from the bonding mode, as shown in 



[figure] [4] []4;b,c). The remaining part of the co+ branch is almost independent of d. This part of 



the excitation branch mimics the single layer plasmon subjected to Eg (symmetric configuration, 
[figure] [4] []4| ^b.4, c.4)) or gap independent (asymmetric case, [figure] [5] []5 (b.4,c.4)). The latter 
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case is of special interest since it obeys the conventional plasmon dispersion law for graphene 
while seemingly progressing into the gap within the particle-hole mode region. However, its spec- 



tral weight as seen in [figure] [7] []7 demonstrates strong Landau-like damping once it enters region 
CO + Vfq > Iji/h. 

By bringing the layers close to each other, CO- becomes acoustic like (see ??), and ,for small 
inter-layer distance, may be mostly accommodated by region 4<. The (0+ mode is similar to that 
for the single layer. That is, (0+ ~ (Op{Eg) for the symmetrical case ?? , and is gap independent in 
the asymmetrical case with co+ ^ 03p{Eg = 0). 



IV. Concluding Remarks 

Theoretical results reported here were carefully presented to stimulate experimental verification 
for both the single and double layer configuration. One possible probe which may be employed 
is EELS using an electron spectrometer In EELS, some of the electrons undergo inelastic scat- 
tering, losing energy and having their paths slightly modified. The inelastic interactions include 
phonon excitations, inter band and intra band particle-hole excitations, plasmon excitations and 
Cerenkov radiation .^^ A separate study is needed to calculate the stopping power due to these 
mechanisms in order to supplement the results reported in this paper. 
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Figure 1: (Color online) Schematic illustration of the pump-probe consisting of circularly polar- 
ized light and a beam of energetic charged particles. On the left is the energy band structure for 
conventional Dirac fermions near the K-point. On the right is the corresponding band structure for 
dressed Dirac fermions. 
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q/Ap q//CF 

Figure 2: (Color online) Labeled regions for analytic results for the polarization function in the 
case of conventional graphene (zero energy gap, left panel) and for the case when there is an 
induced gap Eg / ji = 1.5 (right panel). 
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(a.l) (b.l) 




Figure 3: (Color online) The plasma excitations and their damping are presented for a single 
graphene layer. In rows (1, 2, 3), the energy gaps chosen are Eg //i = {0.0, 1.0, 1.5}, respectively. 
In the panels labeled (a), we plot 3m n('^)(^, co) as well as 9te £{q, (o) = (with d — ^ °°). Panels 
(b) were obtained by plotting the spectral weight — (l/;r)3m Ti{q,(0 — ijq) which involves the 
imaginary part of the RPA form of the polarization function. In panels (a), the darker the shaded 
regions, the larger is the value of 3m Tl^'^\q,0}), corresponding to Landau damping. The black 
solid curves in panels (a) show the plasmon dispersion obtained by solving 9le £{q, (o) = 0. The 
finite lifetime of the plasmon modes represented by the thin solid black curves in panels (b), was 
included in the polarization function to take account of phonon-like scattering and was determined 
by solving 3m Yl{q^ (O + iYq) = 0. The auxiliary dashed lines show n> and n< regions. The upper 
and lower insets of (a.3) correspond to the fixed values of q so that the plasmon branch is below 
and above co = hvpq line. The green dot indicates the position of the plasmon resonance (Op. Red 
(1) curve indicates 3n(°), while blue (2) curve stands for 311. 
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Figure 4: (Color online) Plasma excitations and damping (analogous to column (a) of [fig 



ure] [3][]3| ) for double layer configuration of graphene layers. The imaginary part of the polar- 
ization n^°) is plotted as a function of frequency and wave vector. The solid black lines are the 
plasmon dispersion curves. Columns (a, b,c) are for Eg/ jj. = {0.0, 1.0, 1.5}, respectively. The in- 
duced gap and chemical potential are the same for both layers (symmetric case). Rows (1, 2, 3, 4) 
depict plasmon branches for inter-layer distances kpd = {0.1, 0.5, 1.0, 5.0}, respectively. 
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Figure 5: (Color online) Plasma excitations and damping (analogous to column (a) of [figure] [3] [] 3 



for graphene double layers. The difference with [figure] [4] []4 is that the energy gap Eg may not be 
the same for both layers. Columns (a, b,c) correspond to an induced gap Eg/ ii = {0.0, 1.0, 1.5}, 
respectively, only for one of the layers. Although the chemical potential is the same for both 
layers, a gap is not induced on the second layer (asymmetric case). Rows (1, 2, 3, 4) show plasmon 
branches for inter-layer separations kfd = {0.1, 0.5, 1.0, 5.0}. 
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Figure 6: (Color online) Plots similar to panels (b) in [figure] [3] [] 3 showing 3m n^'j'' for double 
graphene layer configuration. For columns (a,b, c), we chose Egj ji = {0.0, 1.0, 1.5}, respectively. 
The induced gap and chemical potential are the same for both layers (symmetric case). Rows 
(1, 2, 3, 4) display plasmon branches for inter-layer distances kfd = {0.1, 0.5, 1.0, 5.0}. 
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Figure 7: (Color online) Plots similar to those in column (b) of [figure] [3] []3 showing 3m Y^^^ for a 
pair of graphene layers. Columns (a, b, c) were obtained using Egj \x = {0.0, 1 .0, 1.5}, respectively. 
Although the chemical potential is uniform for both layers, a gap is not induced on the second 
layer (asymmetric case). Panels (1, 2, 3, 4) show plasmon branches for inter-layer distances kfd = 
{0.1,0.5, 1.0,5.0}. 
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